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Introduction 

The Possibility of major changes In the global environment pre- 
nts a difficult task to the scientific research community: to devise 

s W hTs°asTh yZing the CaUSCS and Pr ° jectin g the courees of these 
shifts as they are occurring. Purely observational approaches are 

inadequate for providing the needed predictive or anticfpatory infor- 
mabon because response times of many terrestrial ecosystems are 
ow. and there is a great deal of variation from place to place Fur- 
thermore. many important processes, such as soil processes, cannot 
be measured directly over large areas. We need models to express 
our understanding of the complex subsystems of the earth P how 
they interact, and how they respond to and control changes in cli- 
mate and biogeochemical cycles. 6 

The primary research issue in understanding the role of terres- 
nal ecosystems in global change is analyzing the coupling between 
processes with vastly differing rates of change, from photosynthesis 
to community change. Representing this coupling in models is the 

Sh“ enge to modellng the te,Test,,al « pS or £ 

Terrestrial ecosystems participate in climate and in the biogeo- 
chemical cycles on several temporal scales. Examples of processes — 
that operate on short time scales (i.e., less than days) are the meta- 

ce°rtairmfr Se hi r r POnSlble ^ 311(1 mal "tenance, and 

certain microbial processes associated with dead organic matter 

decomposition. The associated energy balance is also affected at 
short time scales. 
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some of the carbon fixed by photosynthesis is incorporated into 
pl^t tissue and is delayed from returning to the atmosphere unttit. 

Is oxidized by decomposition or fire. This slowe ( • ’ y 

months) carbon loop through the terrestrial component of the car- 
bon cvcle which is matched by cycles of nutrients required by 
plants and decomposers, affects the increasing trend m atmosphere 
COo concentration and imposes a seasonal cy 
Moreover this cycle includes key controls over biogenic trace ga 
production. The structure of terrestrial ecosysl Lems which 
on even longer time scales (annual to century), is the integrated 
response to the biogeochemical and environmental constraints that 
develop wer the intermediate time scale. The loop is closed back to 
the climate system since it is the structure of ecosystems, taduding 
{mectesT' composition, that sets the terrestrial boundary condition in 
the climate system through modification of surface roug 
1U j onH tn;* trreat extent, latent heat exchange « 

sca.es contain explicit feedback loops 

wh^h may modify ecosystem dynamics and linkages betw«n ecosys- 
r e mfa^d t”e atmosphere. Consider again the coupling of long-term 
climate change with vegetation change. Climate change win a 
vetretation dynamics, but as the vegetation changes quantify or 
type of structure, this may feed back to the atmosphere by changing 
Ser, energy, and gas exchange. Biogeochemical cycling will also 

change altering the exchange of trace gas species and nutrient aval 
change, auei g . 11m9tp resulting from increased 

ability. The long-term change in climate, resu g 

atmospheric concentrations of greenhouse gases (e.g., C0 2 , CH 4 an 
"oxide INjOl) «m further modify the ^ ^ 

potentially induce further ecosystem change. Modeling these tat 
tions requires coupling successional models to biogeochem c 
^ o Xtalogical models that describe the exchange of water 
energy' and biogenic trace gases between the vegetation and the 
atmosphere at fine time scales. There does not appear to be any obvi- 
ous way to allow direct reciprocal coupling of atmospheric gene 
circulation models (GCMs), which inherently run with fine time 
steps to ecosystem or successional models, which have coarse te 
oSres^tTn, without the interposition of physiologic^ canopy 
models This is equally true for blogeochemlcal models of the 
“xchaLge of carbon dioxide and trace gases. This coupling across 
ume scLs is nontrivial and sets the focus for the modeling strategy. 


Scales of Interactions 

Based on current model structures, atmosphere- Wosphem Inter- 
actions can be captured with simulations operating with three char 


W3>rn<— 3JO o | >s 


Brands Bretfierton et of. 


183 


acteristic time constants (Figure 1). The first level represents rapid 
(seconds-day) biophysical interactions between the climate and the 
biosphere. The dynamics at this level result from changes in water, 
radiation, and wind and accompanying physiological responses of 
organisms. Dynamics at this level occur rapidly relative to plant 
growth and nutrient uptake, and far more rapidly than species 
replacement can occur. Simulations at this level are required to pro- 
vide information to climate models on the exchange of energy, water 
and C0 2 . Tests of this level of model can be accomplished using 
experimental methods including leaf cuvettes, micrometeorological 
observations, and eddy correlation flux measurements. 

The second level captures important biogeochemical interactions 
This level captures weekly to seasonal dynamics of plant phenology 
carbon accumulation, and nutrient uptake and allocation (Figure 1). 
Most existing models at this level use integrative measures of cli- 
mate such as monthly statistics and degree-day sums. Changes in 
soil solution chemistiy and microbial processes can be captured at 
this level for calculation of trace gas fluxes. Primary outputs from 
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Figure 1. Scales of interaction between climatic interactions and ecosystem properties. 
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this level of model are carbon and nutrient fluxes, biomass leaf area 
index (LAI), and canopy height or roughness. This level of model is 
usually tested in field studies with direct measurements of biomass, 
canopy attributes, and nutrient pools or fluxes. ... , 

A third level of model represents annual to decadal changes m 

biomass and soil carbon (net ecosystem 

age) and in ecosystem structure and composition (Figure 1) Inputs 
are statistical distributions of climate variables and calculated 
indices summarizing the effects of climatic conditions on biomass 
accumulation and decomposition. The outputs inclu e ecosys em 
element storage, allocation of carbon and other elementsamong 
sue types, and community composition and structure. Such models 
are currently based either on individual organisms or on species 
correlations with environment. The former are difficult to app y a 
large scales because of computational and data requirements, and 
considerable work will be required to develop large area imp emen- 
tations. This type of model is validated using a combination o 
process studies, as described above. These processes need to be 
integrated and validated in comparative studies to derive annu 
fluxes. The community composition and population dynamics 
aspects of these models are often validated using paleodata. 

are two scales of spaUal resolution. At fine resolution, eh- 
mate model results are used to drive regional-scale ecosystem 
process models, which are rather mechanistic. Much higher spatial 
resolution than currently available is needed for terres ^ 1 
variables, particularly for complex terrain (e.g„ mountain ranges). 
Next-generation GCMs, such as the community climate model ver- 
sion 2 (CCM2) from the National Center for Atmospheric Research 
(NCAR) (with 250 x 250 km resolution), will improve on this p 
lem, but only in a modest fashion, and will not satisfy all needs of 
terrestrial modeling. Development of nested model techniques 
greater computing power will be needed to produce climate-related 
results in the range of resolution (10-20 km) required by ecosystem 
process models. Spatial distribution of precipitation “ 
and improved surface topographic definition (as done in COM2) is a 
critical step in improving orographically related P red P^° n „ ^ 
At the coarser scale required for the operation of global tenresWa^ 
vegetation models, GCM grid scale is nearly adequate, alth g 
improvements are highly desirable. However, any d^amic process 
model of global vegetation can be linked to key satellite^ 
a bies such as advanced very high resolution radiometer (AVHRR) 
vegetation index data and the normalized difference vegetation index 
(NDV1), that provide reference characteristics of the land surface. 
Current global NDVI-based land cover maps have a resolu ion o 
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about 4 x 4 km; however, for many global applications these data 
could be aggregated to near GCM cell scales. It is important that 
future GCMs be able to interactively incorporate satellite data for reg- 
ular redefinition of surface albedo and of vegetation characteristics 
related to evapotranspiration (ET), such as LAI. Additionally, regular 
monitoring of changing land cover and land use will be important. 

Atmospheric GCMs tend to simulate mean atmospheric condi- 
tions. However, many terrestrial processes, such as biome replace- 
ment, are triggered by occasional extreme events. Although these 
events do not have global significance, they are of unsurpassed 
importance regionally. The meteorological conditions that triggered 
the 1988 fires at Yellowstone National Park will have regional conse- 
quences for the next century. Spring frosts, which occur early in the 
growing season; multiple years of successive drought; floods; and 
hurricanes are other examples of extreme meteorological events that 
have potentially significant ecological ramifications. Some lakes in 
central Australia only contain water a few times per century, when 
extreme precipitation events occur. The utility of augmenting GCM 
results with regional climatological statistics, as suggested by F. 
Bretherton in the appendix to this chapter, may be the most reason- 
able way of providing this data. 


Specific Variables Linking the Atmosphere 
with the Terrestrial Ecosystem 

Precipitation 

Simulation of precipitation requires high spatial resolution (10-50 
km), particularly in complex terrain. Minimum event resolution is 
around 3 mm for regional process simulations. At least daily time 
resolution is essential to differentiate between precipitation that is 
intercepted by vegetation canopies and evaporated immediately and 
precipitation that enters the soil rooting zone, with a subsequent 
residence time of days to months. Also, interception and evaporation 
have only physical controls; the dynamics of soil water uptake are 
controlled by plant physiology in concert with physical drivers. 

There are several ways to reduce uncertainties of local precipita- 
tion from coarse-resolution models: 

• Empirical relationships can be developed between mean precipi- 
tation and other large-scale circulation statistics for the region 
defined by a GCM grid cell and precipitation data at individual 
stations within the grid cell (see appendix to this chapter, “A 
Regression of Atmospheric Circulation to Local Weather”). Such 
relationships can then be used to extract subgrid-scale precipi- 



tation variations from GCM simulations. This is the approach 
used in standard numerical weather forecasts. However, it is 
not clear whether such empirical relationships will hold in a 
changing climate. 

. The planned Tropical Rainfall Measuring Mission (TRMM) will 
allow precipitation to be derived from satellite observations (in the 
microwave). Limiting its value, TRMM is planned to be an 
exploratory mission, i.e., of limited lifetime, and it will be focused 
on the tropics. 

• Algorithms are being developed to extract intensities of convective 
precipitation in the tropics from satellite measurements of outgo- 
ing longwave radiation (OLR) at the top of the atmosphere. They 
employ the fact that colder brightness temperature is correlated 
with higher convection and more intense rainfall. The approach 
looks promising, and it should be extended to the rest of the 
globe. High-resolution (15-minute, ==20-km) OLR data are in the 
archives, as are high -resolution precipitation data at river gauge 
stations (U.S. Geological Survey network); these data are useful 
for testing the algorithms. If global algorithms can be developed, 
they offer the possibility of gaining a self-consistent data set for 
expressing subgrid-scale precipitation in terms of larger-scale 
variables such as mean precipitation or OLR. It is not known 
whether the relationships will hold in a changing climate. Looking 
at year-to-year variations in the relationships may provide a clue, 
depending on the quality and availability of data sets for regional 

verification. 

. Mesoscale models nested in GCMs (e.g„ Dickinson et al„ 1989) 
promise predictive capability for local precipitation. A GCM is run 
at coarse resolution, and then a high-resolution mesoscale model 
for a region such as the western United States uses the GCM out- 
put as boundary conditions or driving functions and simulates 
the climate within each grid cell of the region. The usefulness of 
this approach depends, of course, on whether cloud physics and 
precipitation dynamics are properly incorporated in the models. 
For global applications, precipitation with a precision of 10 mm 
may be acceptable. 


Temperature 

Canopy-level temperatures are needed to drive evapotranspiration, 
photosynthesis, and respiration computations, preferably at daily 
time scales. Various growing season definitions, ecosystem phenol- 
ogy, etc., are best defined either by integrated daily temperatures 
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(i.e., growing degree day crop forecasts) or by thresholds, e.g., last 
frost-first frost growing periods. The needed accuracy Is about 1°C. 
Substantial spatial variability within GCM cells is caused by topogra- 
phy (slope, aspect and elevation) and by variability in land surfaces 
such as forest vs. cropland, irrigated land vs. desert, and upland vs! 
wetland. Variability in topography and land cover can be described 
statistically, rather than by geographically explicit treatments. 

Monthly or yearly average temperatures are used to drive models 
of soil processes, such as decomposition and N mineralization. For 
these models, soil temperatures at a depth of 10-20 cm are needed, 
as well as surface or canopy temperatures. Additionally, simplified 
primary production models, once they have been “calibrated" by 
daily canopy models, can be used for general global estimates with 
minimal data requirements. GCM grid cell output of daily maximum 
and minimum temperatures is very useful for defining continental- 
scale vegetation phenology and growing seasons. 

Atmospheric Deposition of Nutrients 

A subset of required precipitation data for terrestrial ecosystems is 
wet and diy deposition of atmospherically transported chemicals, 
including nutrients such as the nitrate, ammonium, sulfate and 
phosphate radicals (N0 3 ~, NH 4 +, S0 4 =, and P0 4 -). Acid rain and air 
pollutant deposition effects on terrestrial ecosystems have been widely 
studied. GCMs simulate the physical mechanism of transport of these 
aerosols, particulate matter, etc., but they do not couple the com- 
pounds to measured source fields, and they do not include any of the 
atmospheric chemistry involved in their transformations. Future ver- 
sions of GCMs should be able to provide deposition estimates globally. 

N0 3 - is deposited from the atmosphere both as nitrate in rainfall 
and directly as nitric acid (HN0 3 ). The fractional input from both 
sources appears to be about equal, as some recent studies of HN0 3 
downward fluxes over grasslands have shown. Dry deposition of 
nitrogen dioxide (N0 2 ) is less important because of the smaller rela- 
tive atmospheric abundance. The atmospheric budget of NO x , which 
consists of nitric oxide (NO) plus N0 2 , is approximately 18.5-89.5 
Tg N/yr, as shown in Table 1. 

In remote areas, the abundance of HN0 3 is typically 0. 1 ppbv or 
0.2 Tg N, with an average atmospheric lifetime against rainout of five 
days. This gives a budget of 50 Tg N/yr (with an uncertainty of at 
least a factor of two). The lightning source may be dependent on lati- 
tude and can be regionally important (Liaw et al., 1990), but will be 
dispersed since the NO it produces must be oxidized to HN0 3 before 
deposition. Most of the large emissions associated with urban pollu- 
tion are observed to decay rapidly (presumably to the nearby areas). 
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Table 1: Global budget for NO x 


Source 

Tg N/yr 

Stratospheric input (from N 2 0) 
Lightning 

Soil sources 
Fossil fuel burning 
Biomass burning 

0.5-1. 5 
8 

(2-20)* 

(81 ± 65.7]** 

8 

(1-16) 

21 

(14-28) 

12 

(1-24) 

Total 

18.5-89.5 

excluding large lightning values 

Logan, 1983; Liu and Cicerone, 1984 

•Liawetal., 1990 

1 


Accurate global modeling of the nitrate system is difficult and has 
only recently been attempted in global chemical transport models 
(CTMs). The current models cannot be used to predict nitrate input 
accurately, but in the next few years they are expected to produce a 
good but incomplete global picture of the NO x budget. Anthro- 
pogenic inputs are expected to dominate near industrial regions. 

Ammonia (NH 3 ) deposition from the atmosphere is most likely to 
come in the form of ammonium nitrate (NH 4 NO 3 ) or some other 
neutralized ammonium aerosol. Although elevated concentrations o^ 
NH 3 gas have been measured above biologically active hot spots 
(e.g., animal feedlots), concentrations are difficult to detect in the 
free troposphere. The implication (consistent with the observed 
abundances of HN0 3 and acidic aerosols) is that NH 3 is removed 
from the atmosphere in aerosols. (The lifetime against oxidation o 
the hydroxyl free radical, OH, is greater than 30 days, even in the 
tropics.) If NH 3 forms aerosols in the boundary layer in the immedi- 
ate vicinity of where it is emitted, then the most likely effect of 
ammonia volatilization is the horizontal dispersion of nitrogen. 

Ammonia releases are deposited within a few hundred to a thou- 
sand km (generally downwind) of their source. Global models or 

atmospheric NH 3 are not available. 

The sources of sulfates deposited to the ecosystem Include marine 
sea salt, photochemically oxidized marine sulfides, and anthro- 
pogenic sulfates. (We exclude here volcanic emissions.) In regions 
where marine sulfide- to- sulfate sources may be important, it is 
likely that deposition of sea salt sulfates provides adequate fertiliza- 
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tion. These regions are thus independent of the atmospheric chem- 
istry. Large perturbations to sulfate-limited ecosystems will proba- 
bly occur only if they are downwind from industrial sources, such as 
areas with considerable combustion of sulfur-containing fuel or 
smelting. 

For ecosystems that are limited by phosphorus, the only atmos- 
pheric source is associated with large dust storms. Such storms are 
not regular, annual processes, but rather extreme events that can- 
not be predicted from current atmospheric CTMs. Changes in inputs 
of phosphorus are likely to be the result of extreme climatic events. 

The potential impact of acid rain is proportional to the total flux 
of hydrogen ions (H+) both in rainfall and in dry deposition (e.g., 
HN0 3 ). Normal rain is acidic: pH’s of about 5.5 are due to the solu- 
bility of 350 ppm of C0 2 in water. Even in remote regions, organic 
acids such as formic acid (the result of oxidation of CH 4 and other 
hydrocarbons) reduce the pH of rainfall to 5 or less. In regions per- 
turbed by large anthropogenic sources of acid precursors (e.g., NO 
and sulfur dioxide, or S0 2 ), the acidity of the rain is likely to be 
below a pH of 4. However, in regions with alkaline soils, the natural 
dust can neutralize the acidity of rain (pH > 6). Regional acid depo- 
sition models have been used to predict H+ fluxes for part of the 
United States, but the deposition of H+ is difflcult to predict on a 
global scale. Significant amounts of acid rain (and deposition) are 
expected a few thousand km downwind of large anthropogenic 
sources of NO and S0 2 . 


Ozone 

The interaction between ozone and vegetation is an important link 
in modeling both the ecosystem and the atmosphere. High local 
ozone concentrations can be deleterious to vegetation, damaging 
stomata and affecting evapotranspiration. Currently there are 
regional models for air quality that predict tropospheric ozone over 
part of the United States. Many research groups are developing 
CTMs for global tropospheric ozone, and we may expect such mod- 
els to contribute to these regional atmospheric-ecosystem models in 
the next few years. The occurrence of extreme ozone events will still 
be difficult to predict in the global models. 

The destruction of ozone at the earth’s surface is an important 
part of the global budget for tropospheric ozone. Observations show 
that ozone is removed with high efficiency above regions with active 
plants. Emissions of NO, N 2 0, CH 4 , isoprene (C 5 H 8 ), and other 
hydrocarbons from vegetation and soils play a major role in global 
atmospheric chemistry. 
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Humidity 

GCMs provide a water mixing ratio in the lowest layer of the 
model (900-1000 mb, depending on the model). In some boun ary- 
laver models, this can be translated into dewpoint at the canopy 
surface Humidity is critical for modeling ET and stomatal response 
models, but less so for soil carbon or nitrogen cycling except as a 
component of hydrologic balance. If a spatially and diumally conser-^ 
vative measure of humidity is used, such as dewpoint or mixing 
ratio fairly broad regional average conditions may be adequate. The 
needed accuracy is - 1°. However, canopy- or surface-level estimates 
are needed, as a bulk tropospheric mixing ratio probably underrep- 
resents the near-surface humidity over vegetated surfaces. 

Solar Radiation . 

Incoming shortwave or photosynthetically active radiation (PAR) is 
very critical for photosynthesis and ET models, although less so for 
modeling soil processes. Incoming solar radiation is well modeled by 
GCMs for a molecular atmosphere, but inadequate simulation o 
clouds and aerosols in current climate models makes this a feature 
of questionable value. Hourly or daily resolution is needed for photo- 
synthesis and ET models; monthly or annual totals are used n 
ecosystem dynamics models. It Is not critical to partition between 
direct and diffuse radiation or between total radiation and PAR. 
Subgrid-scale topographic variability can easily be handled with cor- 
rections for slope and aspect. However, subgrid-scale clouds and 
optical variability are problems. 

Wind and Dust 

Extreme storms, hurricanes, etc., are important ecosystem s- 
turbance triggers. The recurrence of extreme climatic events deter- 
mines ecosystem structure to a certain extent. There rea y s no 
satisfactory way to simulate extreme winds with a coarse-resolution 
model. Parameters such as minimal surface pressure may e use 
ful for defining storm existence and perhaps storm tracks. How- 
ever, an operational diagnostic of maximum winds may not be 
readily available in GCMs. Average wind velocity is also an impor- 
tant component of evaporation, and daily averages may be ade- 
quate for ET estimates. This parameter, average wind velocity s 
highly variable and unpredictable in complex topography t a s 

below grid-scale resolution. , x . 

Extreme wind conditions are of high importance for dust plumes 
and soil erosion calculations (which also require soil texture and 
moisture information). Dust affects soil properties, incident solar 
radiation, and terrestrial productivity. All of the necessary atmos 
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pheric parameters are available in principle in GCMs. However, the 
current capabilities of the GCMs and our understanding of source 
functions need significant improvements before useful applications 
can be made. 

Snowpack 

Snowpack is a critical available output from GCMs. Snowpack is 
a major albedo and hydrologic balance determinant and is an 
important output variable of most GCMs. Obviously there is much 
subgrid-scale, topographically induced variability in snowpack 
dynamics. However, even grid-scale snowpack information, as pro- 
vided currently by GCMs, is useful. Subgrid-scale modeling is pos- 
sible from higher-resolution topographic and topoclimatological 
models, and the U.S. Soil Conservation Service operates a national 
Snow Survey Network (used for summer irrigation scheduling) of 
ground measurements for model development /validation. 

Snowpack is a critical mechanism for water storage, providing 
summer water in arid regions and controlling flood potential in 
many wetter areas. Snowmelt is an important trigger of growing sea- 
son, causing an intense burst of trace gases from the spring soil and 
hydrologic flushing of dissolved elements. Spring snowmelt initial- 
izes calculations of seasonal soil water depletion and ecosystem 
stress and produces primary annual hydrologic discharge in semi- 
arid lands. Snowpack insulates soils from extreme temperatures, 
and is important in vegetation/crop survival. 


Variables Linking Terrestrial 
Ecosystems to Atmospheric GCMs 

Evapotranspiration 

The single most important feedback from the land to the climate 
models is the partitioning of incident solar radiation (H) to sensible 
or latent heat (LE), the Bowen ratio (H/LE). This partitioning is con- 
trolled by the surface evapotranspiration in a complex fashion and 
is computed with soil-vegetation-atmosphere models. Major 
decreases in seasonal Bowen ratio occur as vegetation develops, 
increasing ET and latent energy. Later in the growing season, as 
vegetation either senesces or endures soil water deficits, physiologi- 
cal water stress and stomatal closure can cause Bowen ratios to 
increase by an order of magnitude as progressively more of the 
incoming energy produces sensible heat. This cannot be effectively 
simulated by meteorological measures alone, but can be calculated 
(interactively) in models of ecosystem water balance that determine 
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seasonal LAI, soil moisture, and water use efficiency as a function of 
climate C0 2 concentrations, plant physiology, and other perturba- 
tions and stresses. Ecosystem water balances are typically executed 
on hourly or daily time scales. 


At minimum, ecosystem models predict seasonal phenology and 
assume some ecosystem structure and soils. A look-up table trans- 
lation scheme can be used to yield broadband albedos for each of 
these ecosystem compartments, which can then be weighted o g ve 
the mean albedo for a GCM grid cell as a function of time. There are 
sophisticated models (e.g., ray tracing models) that calculate spec- 
tral reflectivities as a function of, among other things, leaf shape 
orientation, and distribution. However, significant generalization of 
these models is needed before they are useful for GCM-scale appli- 
cations. Inversion of satellite -derived surface reflectance can now 
provide global maps of surface albedo. 

A major uncertainty in GCMs is the highly dynamic albedo 
changes caused by the delivery and melting of snow. Albedo varies 
with age (and depth) of snow and with the masking depth of vegeta- 
tion. This introduces a large degree of uncertainty to the energy bal- 
ance of high latitudes. Either mechanistic submodels of snowpack 
dynamics or regular reparameterizations of the GCM during simula- 
tions by satellite data are probably required to meaningfully improve 
this situation. 


Roughness Length 

Topographic roughness and vegetation height affect momentum 
dissipation at the surface. In general, topographic roughness domi- 
nates the momentum exchange. The improved spatial resolution o 
the NCAR CCM2 allows higher definition of grid-cell topograp y. 
However, additional improvement may be possible by using subgrid- 
scale information on topographic variation. 

Typically, a roughness length is assigned to each vegetation type, 
ranging from around 0.02 to 6.0 m, small relative to topography. 
Models that simulate ecosystem transitions can, in principle, output 
roughness length as a function of time. How important this is for 
altering GCM circulation has not been established. 


Greenhouse Gases 

Models of canopy photosynthesis-respiration balances an soi 
carbon cycle dynamics in ecosystems by definition keep track of 
surface C0 2 uptake from the atmosphere by simulating the produc- 
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tion and return of C0 2 to the atmosphere through soil organic mat- 
ter turnover. When applied to scenarios of climate change, these 
models project changes in the net fluxes of CO a that are useful for 
studies of the carbon cycle. A daily or longer surface C0 2 balance 
could in principle be entered into a GCM. Again, the prtmaiy prob- 
lem is the spatial aggregation of the ecosystem model outputs to the 
GCM cell size. Alternatively, simple satellite-driven AVHRR/NDVI- 
based models have been highly correlated with C0 2 balance and 
could be implemented globally (Running and Nemani, 1988). 

Carbon dynamics models can be adapted to the study of other 
trace gases. We envisage that a separate treatment would be 
required for each gas. N s O and NO fluxes require elaboration in soil 
compartments and explicit treatment of soil moisture regimes. Iso- 
prenes, terpenes, and other nonmethane hydrocarbons require 
some scaling to the metabolic function of certain vegetation species 
Modeling of carbon monoxide must be linked to a fire model, and 
CH 4 requires a detailed hydrology model. Because of the many link- 
ages required, global-scale modeling of trace gas fluxes is not cur- 
rently possible. 

Unmodeled Land Parameters Needed by GCMs 

A number of soil and vegetation characteristics are needed for 
GCM surface parameterization, as exemplified by the biosphere- 
atmosphere transfer scheme (BATS; see Tables 2 and 3 in Dickinson 
et al., 1986). Many of these characteristics are not outputs of 
ecosystem models but are surface parameterizations also needed by 
them. Some of these parameters can be estimated globally, but 
many are impossible to do accurately. Satellite-derived estimates of 
vegetation cover and seasonal variation are available from weekly 
composite global AVHRR/NDVI maps. These can be translated rea- 
sonably well to maximum-minimum LAI for different biome types. 
Physiological parameters like minimum stomatal resistance or leaf 
light sensitivity cannot be directly inferred. Global albedo is now 
being monitored weekly by satellite and could be entered into the 
GCMs. Soil physical and hydraulic parameters cannot be monitored 
by satellite, and existing global soils maps are taxonomically based 
and of questionable accuracy. It is not clear how this problem can 
be improved (see Parton et al. report, this volume). 
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Appendix: A Regression of 
Atmospheric Circulation to Local Weather 

A critical gap in connecting atmospheric climate models to ecosys- 
tem models of the land surface lies in the provision of local estimates 
of surface precipitation and daily weather statistics. Atmospheric cli- 
mate models are designed to operate on spatial resolutions of several 
hundred kilometers horizontally and 1000-300 millibars vertically 
whereas ecosystem models are developed and tested on data relevant 
to a small watershed, typically of 1-10 km. The variables predic ed 
explicitly in GCMs (pressure, temperature, moisture, and win in e 
free atmosphere) relate to the atmospheric circulation on regional 
scales. The lowest atmospheric cell may be between 1000 and 950 
mb and will not be representative of surface conditions wit ou a 
boundary-layer submodel. Subgrid-scale topography, land cover, and 
lake effects can all have a major impact on local temperature and 
precipitation, particularly the daily maximum and minimum temper- 
atures. the incidence of solar radiation, and the probabilities of rai 
falls of different intensities, which are the variables needed to run 
ecosystem models. Although for biogeochemical purposes monthly 
averages are adequate, to simulate the hydrologic balance, daily data 
and a knowledge of the diurnal cycle are imperative. Successional 
models are sensitive not only to seasonal averages but also to 1S U1 " 
bances such as unseasonal frosts and fires. Giving a general defini- 
tion of disturbance is difficult, but the meteorological factors involved 
can all be derived from a good daily weather record. Though surface 
temperature, solar radiation, and rainfall are all computed hourly in 
a GCM as part of the treatment of the atmospheric boundary layer, 
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these are highly parameterized estimates on a scale of the GCM grid, 
and it is not obvious how to derive locally valid values, 

A similar problem is encountered in weather prediction, in which 
model simulations of the atmospheric circulation have to be trans- 
lated on a daily basis into specific local forecasts for transmission to 
the public. It is believed that changes in local weather are largely 
controlled, at least in a statistical sense, by the regional atmos- 
pheric circulation, so a widely used technique is to correlate over 
several years the values of a selected subset of variables in the 
weather prediction model with the daily record from a local weather 
station. The correlation table is then used, in conjunction with 
future values from the prediction model, to predict the local sta- 
tion’s weather. This technique, known as model output statistics, 
has proved as skillful as most human forecasters using more sub- 
jective methods. 

This approach could easily be adapted to translate the output of 
climate models into equivalent local climates. Running a climate 
model in forecast mode for many years is normally not practical; 
instead, the observed values of the regional-scale circulation vari- 
ables required for the correlation may be taken from actual observa- 
tions rather than from a forecast, using a daily weather analysis. Of 
course, it is not to be expected that subsequent application of this 
correlation table to the output of a climate GCM will automatically 
produce a good simulation of local climate if the GCM produces 
inadequate statistics of the regional circulation patterns. However, 
the exercise will certainly focus attention on which aspects of the 
GCM circulation are Incorrect, and, in conjunction with a local 
ecosystem model, could lead to a figure of merit by which Improve- 
ments in simulations of present climate and the significance of pro- 
jected changes could be assessed. 

An issue that at once arises in applying this approach is the 
availability of suitable data. In most regions of the world there are 
many weather stations recording daily maximum and minimum 
temperature, total rainfall, and hours of sunshine. However, long 
records of hourly precipitation are less readily available, and inten- 
sity of solar radiation (as opposed to the more qualitative surrogates 
of cloudiness and hours of sunshine) is measured at only a few sta- 
tions. However, surface solar radiation can be inferred fairly directly 
from geostationaiy satellite measurements of solar reflectance. Cor- 
rections must be made for variations in viewing angle, for surface 
albedo, and for absorption of solar radiation within clouds and other 
aerosols, and the results must be calibrated against the available 
surface radiation measurements. Once this is done, the data 
required are widely available over most, though not all, of the world. 
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Another Issue is whether the structural relations between a given 
regional circulation pattern and local weather are likely to be stable 
under conditions of changing climate. To the extent that these con- 
nections are controlled by the interactions of topography, wind 
direction, and temperature through the troposphere, changes are 
indeed likely to be caused by changes in frequency of a regional- 
scale pattern and allowed for in this methodology. To the extent that 
they depend upon changes in aerosol or cloud condensation nuclei, 
which are not discussed in routine weather analyses, the correlation 
tables would change over decades. Indeed, systematic trends 
observed in such correlation tables might be the best indicators of 
this type of climate change. 

It is also clearly inadequate to be able to estimate local climates only 
at points where there are weather station records. Modeling of the 
effects of topography and other causes specific to the local weather sta- 
tion is clearly an integral part of any application. This modeling may be 
empirical or may be based on more sophisticated local process studies 
using dynamical mesoscale models for phenomena such as lake effects 
and mountain drainage. In some locations such studies are already 
going on for reasons other than climate change, but it is important to 
generalize this experience to a global scale. 

There are also extreme events of great importance to certain 
ecosystems, such as hurricanes and tornados, of which the fre- 
quency may well change with climate regime. Yet these events are 
not predicted within present climate models. Because of their rarity, 
such events are unlikely to appear in daily correlation tables, 
though special studies might well reveal the regional circumstances 
under which they are most likely to occur. GCM simulations should 
then be examined for changes in hurricane and tornado potential, 
which would be fed offline into ecosystem models. 

If applied to a suitable sample of weather stations within a region, 
the correlation technique would also provide data sets to test and 
validate the physically based parameterizations of the surface vari- 
ables predicted in the climate model for the purpose of keeping 
track of the surface energy and moisture balance. Indirectly, these 
parameterizations depend on the same regional atmospheric vari- 
ables. It is, of course, straightforward to produce similar correlation 
tables within the climate model itself. It could be argued that it is 
inconsistent to maintain within the same framework two distinct 
connections between regional circulation variables and surface cli- 
matology and water balance, one empirical and site- specific, the 
other physically based but highly parameterized. In the long run 
this argument is probably correct, but until more experience has 
been gained the dualism is likely to be a source of strength. 
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